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Abstract 

We construct explicit integral relations between propagators of generalized Schrodinger equa- 
tions that are linked by higher-order supersymmetry. Our results complement and extend 
the findings obtained in [5] for the conventional Schrodinger equation. 
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1 Introduction 

This note is directly motivated by the recent work [9] , where surprisingly simple expressions for 
propagators of supersymmetry-related Schrodinger equations were constructed. In general, the 
quantum-mechanical supersymmetry (SUSY) formalism connects two Schrodinger equations for 
(usually) different potentials, such that solutions of one equation are mapped onto solutions of 
the second equation [3]. The mapping used for relating the solutions to each other is called 
SUSY- or Darboux transformation. While the Darboux transformation has been around for 
more than a hundred years [4j, its application to quantum- mechanical problems within the 
SUSY context is much younger, starting with its relation to the Infeld-Hull factorization [6]. 
Ever since then, the quantum-mechanical SUSY formalism has been extensively used to study 
solvability and spectral properties of Schrodinger equations, for an overview the reader may 
consult [3] [7] [15]. Besides the fact that Schrodinger equations related by the SUSY formalism 
allow for mapping their respective solutions onto each other, such equations are connected in 
more ways. In particular, their propagators and the traces of their Green's functions are linked 
by relatively simple formulas, as was shown in [9] and [TO] [13], respectively. It turns out 
that the link between the Green's functions even persists under linear generalizations of the 
Schrodinger equation [12], such as the effective mass case or minimal coupling to a magnetic 
field. Motivated by this result, in the present note we study the question whether the propagator 
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relation found in [9] can be extended to generalized Schrodinger equations. While for first-order 
SUSY transformations a positve answer has been given [11], here we will study higher-order 
transformations. Section 2 gives a short review on the necessary preliminaries, while in section 
3 we construct our propagator relations for generalized Schrodinger equations. 



2 Preliminaries 

In the following we briefly summarize basic facts about generalized Schrodinger equations, the 
SUSY formalism, propagators and Green's functions. 



The generalized Schrodinger equation. We consider the following generalized Sturm- 
Liouville problem on the real interval (a, 6), equipped with Dirichlet boundary conditions: 

fix) ip"(x) + f'{x) rp'(x) + [E h(x) - V{x)\ ij)(x) = 0, x€(o,6) (1) 

V(a) = m = 0. (2) 

Here /, h, V are smooth, real functions, with /, h positive and bounded in a and b. The 
constant E will be referred to as energy, and in solutions of (P), (|2|) that belong to the discrete 
spectrum, E stands for the spectral value. Any solution i\) of ([I]), ([2]) belonging to a value E from 
the discrete spectrum, is located in the weighted Hilbert space L|(a,6) with weight function h 
[5]. The lowest value of the discrete spectrum will be called the ground state and denoted by 
Eq with corresponding solution ipo. The interval (a, b) can be unbounded, that is, a or b can 
represent minus infinity or infinity, respectively (however, if a and/or b are finite, then we require 
/, h, V to be continuous there). We see that the problem (pQ), ([2]) can be singular, which means 
that its spectrum can admit a continuous part. Equation ([1]) will be referred to as generalized 
Schrodinger equation, since its special cases are frequently encountered in Quantum Mechanics, 
such as the Schrodinger equation for effective mass or with a linearly energy-dependent potential. 
In the quantum- mechanical context, E denotes the energy associated with a solution t/j, and V 
stands for the potential. 



Generalized SUSY formalism. We will now summarize basic facts from p3], for details see 
the latter reference. The boundary- value problem ([I]), ([2]) can be linked to another problem of 
the same kind by means of the SUSY transformation method. Consider 

f{x) V'ix) + f\x) W'ix) + [E hix) - Uix)] *(a?) = 0, x G (a, b) (3) 

*(o) = ¥(&) = 0, (4) 

where the same settings imposed for (jTJ, (|2J) apply. Clearly, a solution ^ = ^(x) and the 
potential U = U (x) are in general different from their respective counterparts ip and V. Now, 
suppose that ip and Uo, ...,u n _i are solutions of the boundary-value problem ([1]), ([2]) and of 
equation ([1]) at real energies E and Xj < E, j = 0, ...,n — 1, respectively. Define the n-th order 
SUSY transformation of tjj as 



m 

hix) 



Wju ,...,u n - 1 ,'ilj)ix) 
W(«o,-,u„_i)(x) ' 1 J 



where W denotes the Wronskian of the functions in its argument and the upper index x of 
D denotes the variable which the derivatives in the Wronskians are applied to. The function 
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^ = D^ Un _ 1 ip as defined in ([5]) solves the boundary- value problem ([3]), 
U is given in terms of its counterpart V as follows: 



if the potential 



U 



+ n /i 
n 2 

+ T 



f f (f 



r\2 



f'h' , 3 / (/» 



A2 



A '/ 
/fc" 



2 / h 2 f h 2h 2 

3 / (h'f 



W? + /' /»' 



2 /i 3 /i 2 



+ 



2/ 



2 /i 2 



h 



(6) 



note that for the sake of brevity the arguments were left out. In the special case n = 1, the 
transformation ([5]) simplifies to 



D x uo Hx) 



70) iy(n ,V)(x) 



«o0) 



TO 
/i(x) 



U (x) 



(7) 



It is well-known [lj [14] that the n-th order transformation ([5]) can always be written as an 
iteration (or chain) of n first-order transformations ([7]), that is, 



D 



uo,...,Un-\ 



n-1 

n 

3=0 



(8) 



where each function Vj solves the equation ([T]) that is obtained after the j-th first-order SUSY 
transformation. Note that (J5]) and ([7]) remain valid when multiplied by a constant, which can 
be used for normalization. Now, depending on the choice of the auxiliary solution u in (J7J), the 
discrete spectrum of problem ([3]), (jlj) can be affected in three possible ways: if A = Eq and 
u = ipo, then Eq is removed from the spectrum of ([3J) , (|3J). The opposite case, creation of a 
new spectral value A < Eq, happens if the auxiliary solution u does not fulfill the boundary 
conditions dU). Finally, the spectra of both problems (UJ, ([2]) and ([3]), (HJ are the same, if we 
pick A < Eq and an u that fulfills only one of the boundary conditions ([2]). 



Propagator and Green's function. The propagator governs a quantum system's time evo- 
lution. For a stationary Schrodinger equation, the propagator K has the defining property 



0) 



exp(-i E t) ip(x) = J K(x,y,t) ifj(y) dy. 

(a,6) 



Suppose problem (UJ, ([2]) admits a complete set of eigenfunctions (ip n ), n = 0, 1,2, ...,M E No, 
where M can stand for infinity, and (4>k), k € R, belonging to the discrete and the continuous 
part of the spectrum, respectively. Then the propagator K has the representation 



K(x,y,t) = h(y) 



M 

y^nQc) exp(-i E n t) ip n (y) + / (f)k(x) exp(-i k 2 t) cj) k (y) dk 

n=0 m 



(10) 



where E n and /c 2 stand for the spectral values belonging to the discrete and continuous spectrum, 
respectively. The Green's function G of the problem (UJ, ([2]) has two equivalent representations 
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[5], both of which we will use here. In order to state the first representation, let tfto i and ijjo,r 
be solutions of equation ([!]) that fulfill the following unilateral boundary conditions: 







0. 



(11) 



The Wronskian W, 



i/>o w>Po r °1 these funtions is given by 

Co 



W^.^ Jx) 



(12) 



where cq is a constant that depends on the explicit form of t/jn i and tpo r . Now we can give the 
first representation of the Green's function Gq for our boundary value problem (pTJ) , $2§: 



G(x,y) 



1 

co 



4>o,i(y) ipo,r( x ) °( x -y) + *Po,i( x ) ^oAv) % - x ) 



(13) 



where cq is the constant from (|12p and stands for the Heaviside distribution. The second 
representation of the Green's function G can be obtained as follows, provided problem (fl]), J2]) 
admits a complete set of solutions: 



M 



G(x,y) 



n=0 



E n -E 



k 2 -E 



(14) 



where the notation is the same as in (jlOp . Note that the Green's function is taken at energy E. 

Propagators related by first-order SUSY. In order to obtain a relation between the 
propagators of the two boundary- value problems ([I]), ([2]) and ([3]), (j4|), we take the propagator 
K\ of the second problem and express it through quantities related to the first problem. For 
the sake of simplicity we assume for now that the two boundary-value problems have the same 
discrete spectrum and that both of them admit a complete set of solutions belonging to a 
discrete and a continuous part of the spectrum. Furthermore, we assume that the solutions of 
problem ([T]), ([5} are real- valued functions. This is no restriction, as equation ([T]) involves only 
real functions. We then find the following relations |llj between the propagators K\ and Kq of 
our boundary- value problems ((TJ) , (EI) , and ([3]), (Jlj), respectively: 



K x {x,y,t) = h(y)D x U0 Dl f K (x,z,t) G (z,y) dz 



(o,6) 



Ki(x,y,t) = h(y) < 



D x D y 



K (x,z,t) G (z,y) dz 



(a,b) 



ifi(x,y,i) = h(y)Dl Dl [ K (x,z,t) lim 



Ga(z,y) 



+ 0_i(ac) exp(-i A t) <f>-i(y) 
^o(z) V'o(y) 



En — E 



dz. 



(15) 



(o,6) 



The first of these relations is valid, if both boundary- value problems admit the same discrete 
spectrum. The second relation applies, if ©, (j4|) admits an additional discrete spectral value A 
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with corresponding solution <f>-\. Finally, the third relation holds, if the initial problem (JTJ) , ([2]) 
has one discrete spectral value more than its transformed counterpart. In this last case, if we 
choose the auxiliary function to be the ground state tpQ of our initial problem, we can simplify 
([15l) as follows: 



Ki(x,y,t) 



'%) 1 



f(y) i>o(y) 



D% 



, J K (x,z,t) V>o0) dz. 
(y,b) 



(16) 



It is immediate to see that the above propagator relations simplify to their well-known conven- 
tional forms [9], if we set f = h = 1. 

3 Propagators related by higher-order SUSY 

From now on we will assume that a higher-order SUSY transformation was applied to the 
initial boundary- value problem flU, fl2J), giving the associated problem (jSJ), Q. We are looking for 
an explicit relation between the propagators of these two problems. To this end, we distinguish 
whether the SUSY transformation adds new values to the discrete spectrum, removes some from 
it or whether the spectra of both boundary-value problems stay the same. 

3.1 Creation of spectral values 

Let us first assume that our iV-th order SUSY transformation creates a N new discrete spectral 
values in the corresponding transformed boundary-value problem. We denote these values and 
their corresponding solutions by E_ n and ^_ n , n = 1,...,N, respectively. According to (fT5j) . 
the propagator of problem (EJ), (|4]) is then given by 



K N (x,y,t) 



Kv) 



M 

^2^ n {x) exp(-z E n t) * n (y) + / <f> k (x) exp(-« k 2 t) dk + 

n=0 ra 



N 



+ V-n(x) exp(-t E- n t) V-niy) 



n=l 



(17) 



Next, we take into account that all functions ^f n , n = 0, ...,M, and k € R, have been 
obtained from solutions by means of an iV-th order SUSY transformation ([5]), using the 
auxiliary solutions Uj, j = 0, ...,N — 1: 



K N (x,y,t) = huj) IK A i Di A , 

+ / L\ 4> k {x) exp(-i k 2 t) 4> k (y) dk 



M 



7,L% ^n{x) exp(-i E n t) tp n (y) + 

n=0 

N 

+ h(y) ^2^^ n (x) exp(-i E_ n t) *_„(y), 



n=l 



(18) 
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where normalization constants and La, were introduced. Before we determine these constants, 
we make use of the defining property ©, transforming (|18p into 



D y 

u ,...,u N _ 1 J ^U ,...,U N _ 1 



K N (x,y,t) = h(y) D, 
+ L l K (x,z,t) <t> k (z) dz (p k (y) dk 



M . 

E L V> I K (x,z,t) ip n (z) dz ip n (y) + 
n=0 (0,6) 



(a,b) 



N 



Hv) ^Z^-n{x) exp(-« E_ n t) *_ n (y), 



n=l 



(19) 



Our normalization constants Lw, and must be chosen as follows [2] 

iV-l / - JV-1 / " 

p=o V p p=o v p 

substitution of which renders our propagator relation (|19p in the following form: 

M N-l 

K N {x,y,t) = h{y)D° m D« m ' 



K Q (x,z,t) ip n (z) dz ip n (y) + 



(o,6) 



Af-1 



H p=0 P (a,b) 

N 

+ h{y) Y^^-nix) exp(-z E_ n t) ^_ n (y), 

n=l 

Next, we rewrite our normalization constants according to the decomposition 

JV-l ^ N-1N-1 ^ ^ 

n = e n X_ - X„ R _ 



(20) 



P =o ^ " X p 



p =0 7=o \~ A P E ~\ 



AT-1 N-l N-l ^ ^ 

II iirr = e n 



p=0 ^ 



We plug this into our propagator relation (|20p and obtain after regrouping terms 
K N (x,y,t) = 

N-l N-l T M 



.a ^o,.,^-, E II xhc \ K °^ z ^ E 

N 



lpn(z) tpn{y) 

E — \„ 



dz + 



+ 



J K (x,z,t) 



<f>k{z) <f>k{y) 
k 2 — A„ 



(a,6) 



+ h(y) ^tf_ n (x) exp(-i£_ n i) *_„(y). (21) 



71=1 



In the final step we make use of the representation (I14p of our Green's function, which converts 
expression (|2ip to the form 



K N (x,y,t) = 

= h(v) D x D y 



N-l 



( X ■ i I r 

E Ht^IT K (x,z,t) G (z,y) dz 
g=0 A i P I J 



+ 



N 



+ h(y) Y,^-n(x) exp(-i E_ n t) <f>- n (y). (22) 

n=l 

Note that the Green's function depends on the index n, as it must be taken at energy A n . 
Expression ([22]) reduces correctly to the conventional result [S], if h = 1 is substituted. 

3.2 Annihilation of spectral values 

We will now assume that our SUSY transformation removes n values from the discrete spectrum. 
Speaking in terms of the factored transformation, each step will remove the corresponding ground 
state from the current system. In other words, after a SUSY transformation of order n, the 
first n discrete spectral values of the initial boundary- value problem will have been removed 
from the discrete spectrum. This particular ordering of the spectral values and solutions to 
be deleted does not constitute a restriction [9], but facilitates notation and calculation. Before 
we turn to the propagator relations, it is necessary to set up some notation in order to obtain 
representations for the auxiliary solutions. In order to do so, we will consider our higher-order 
SUSY transformation (J5]) in its factored form ([8]), that is, an n-th order transformation is 
seen as the n-fold application of a first-order transformation. Each first-order transformation 
yields a new boundary- value problem of the type ([3]), ([4]), such that in total a sequence of n 
boundary-value problems is generated. We will assume that in each transformation step the 
lowest value is deleted from the discrete spectrum, such that after n iterations of our first-order 
transformation the lowest n spectral values are removed. The auxiliary solutions used in each 
iteration of our SUSY transformation will be named Uj where the first index denotes the 
number of the boundary- value problem associated with Uj t k, starting from j = 0. The second 
index in Uj± stands for the solution number associated with the fc-th discrete spectral value. In 
particular, u n ^ n is the ground state of the n-th boundary-value problem and the functions u n j 
for j = n + 1, M represent solutions of the n-th boundary- value problem. Finally, the u n j for 
j = 0, ...,n — 1 stand for solutions of the equation associated with our boundary- value problem 
at energies that do not belong to the discrete spectrum, which implies that they do not fulfill 
the boundary conditions ((4|). For the sake of convenience, let us now introduce abbreviations 
for the SUSY transformation operators ©, using the notation for our auxiliary functions. For 
natural numbers n and j < n — 1 we define 

L n ,j ~ ^u^j,u j>j+1 ,...,u^ n -i- (23) 

This operator maps solutions of the j-ih boundary-value problem onto solutions of the n-th 
boundary-value problem. Therefore, the operator (|23|) admits the contraction property 

L n,j L tk = L n,k- (24) 
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In terms of these operators, our auxiliary solutions u n j for j = n, M are built by the following 
rule: 



u n ,j(x) — u n -i t j(x). 



(25) 



The remaining auxiliary solutions u n j for j = 0, ...,n — 1, which do not fulfill the boundary 
conditions, are constructed as follows: 



u j(x) 



L n,0 U ,j(x) 
Uqj (x) 



f(x) ul ^x) 



dx. 



Note that the functions uoj and uqj solve the same equation and are linearly independent [8]. 
Now we are in position to construct closed-form expressions for our auxiliary solutions. To 
this end, we take the well-known expressions from the conventional case / = h = 1, which 
were obtained in [9], and map them to the present, generalized context. Let v n j stand for the 
conventional auxiliary functions, where the indices have the same meaning and vary exactly as 
in our generalized case described above, then we have from [9] 



Vn,j{x) 



n-1 



Ej) (E n _2 — Ej)...(Ej + i — Ej) 



Wnjvpfi, ...,V 0t n-i)(x) 
W(v 0fi , ...,V , n -l)(x) ' 



(26) 



where the modified Wronskian W n is obtained from W by removing the (n + l)-th row and 
column from the underlying matrix. We will now rewrite the auxiliary solutions v n j and their 
Wronskians in (|26p in terms of the present auxiliary functions u n j. In order to do so, we make 
use of the following results taken from [14] : 



v n ,j(x) 

W n (v fi, ...,V , n -l)(x) 

W(v 0fi , ...,v n ^ 1 )(x) 



[f(x) h(x)}* U n j(x) 

(n-1) 2 (n-1) (n-3) 
[f(x)}^ [h( X )\ i W n ( )(x) 

n(n — 2) 

[f(x)}— [h[x)\ — W{uo,o,-,u ,n-i)(x)- 



We plug these equalities into (|26j) . which then becomes 

i 

[f(x) h(x)Y u nJ (x) = 



(E n -i — Ej) (E n _2 — Ej)...(Ej + \ — Ej) 



1 



h(x) 



h(x) 

m 



n_l 

2 4 W n (u 0fi , ...,no in _i)(x) 

W(u 0fi , ...,U , n -l)(.x) 



Using the abbreviation 



a 



{E n -i — Ej) (E n _2 — Ej)...(Ej + \ — Ej) 



renders our auxiliary solution u n j in the following form: 



h(x) 



h(x) 
W) 



2 W n (u 0fi , ...,U 0) n-i)(x) 

W(u ,o, ...,U , n -l)(x) ' 



(27) 
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We will now show that the propagator K n of the boundary- value ([3]), ([ID, which is obtained 
after an n-chain of SUSY transformations, can be given in the followin form: 



K n (x,y,t) 



Kv) 
M 



n-1 



■^^oESS / K (x,q,t) uo tj (q) ,hi. (28) 



3=0 



W{y) 



(y,b) 



where Kq is the propagator of the initial boundary- value problem ([T]), ([2]). In order to establish 
relation (j28f) . we will use induction, proceeding similarly to how it was done in the conventional 
case [9]. The first-order case K\ has already been established [TT] and is given in ([TO]) . Now 
assume this relation to hold for K n , then our induction step starts at 



K n+1 (x,y,t) 



'%) 1 



f(y) u n>n (y) 



Ln+l t n J K n (x,Z,t) U n>n (z) dz. 



(29) 



(y,b) 



Since we assume that ([25]) is true for n, we can substitute it into ([29"]) . After ordering terms, we 
get 



K n+l (x,y,t) = 
V /(y) "n,n(y) 



n-1 



3=0 



h(z)\v Wj(z) 



(y,b) 



f(z)J W(z) 



u n ,n{z) J K (x,q,t) u j(q) dq dz. 
(z,b) 



Next, we substitute the ratio of Wronskians Wj/W by means of our representation (|26|) . make 
use of the contraction (|24p and arrive after some simplification at 

K n+1 (x,y,t) = 



f(y) u n ,n{y) 



L n+ifl^2~F, / K z ) u n ,n(z) u ntj (z) / K (x, q,t) u ,j{q) dq dz. 

1-° (y,b) (z,b) 



The area of integration forms a triangle T in z-g-space, which can be seen as the upper half of 
the rectangle [y, b] x [y, b]. Therefore, integration over T can be replaced by integration over the 
rectangle minus integration over its lower triangle: 



K n+1 (x,y,t) 



'%) (-1) 



n-1 



n-1 



f(y) u n! n(y) 



Ln+1,0 



(-1) 3 ' 

,=o Cn >i 



h(z) u n , n (z) u n j(z) dz J K (x,q,t) u j(q) 
(y,b) (y,b) 



K (x,q,t) u j(q) J h{z) u n>n (z) u n j(z) dz dq 

{V,b) (9,6) 



(30) 
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Before we process this expression further, let us evaluate the integral with respect to z: 



h (z) Un , n ( Z ) u nJ (z) dz = m w ^ z ^m) _ m wp^m (31) 



Ej — E n 



Ej — E n 



This can be verified in a straightforward way by diffentiating the right hand side and replacing 
the second derivatives by means of our equation ([1]). The second term on the right-hand side 
of (]3ip is a constant and will cancel out in subsequent calculations, whereas the first term will 
now be modified further. First we make use of the relation 



replacing the Wronskian in (|3 1 [) : 



7(0 W{u n ^U n M) 
HO U n ,n(Q 



h{z) U ntn (z) U n j(z) dz 



V/(0 MO «».n(0 4w,„«nj(0 f(P) W{u„ n ,u n ;){b) 



Ej — E n 



Ej — E n 



Ej — E n 



Ej — E n 



(32) 



We now replace the function u n +ij by its representation (|27p for n -\- 1 and simplify the result: 



h{z) u n , n (z) u n j(z) dz 



MO 
/(0 



)(0 f(b) W(u n , n Kb) 



W(u ,o, ...,UQ,n)(0 



Ej — E n 



(33) 



We can replace the integral (|33p in our propagator relation (|30p . note that the integral ap- 
pears twice there. After some simplification and regrouping terms, we arrive at the following 
expression: 



K n+1 (x,y,t) 



n+l 



l) n L* n+li0 J2(-iy J K (x,q,t) Uo,M dq + 

■? =0 <„ M 



+ 



'%) ("1 



\n— 1 



f{y) u n ,n(y) 



En+1,0 



K (x,q,t) 



%) 

m 



(y,b) 

u n ,n(q) 
W(q) 



n-l 



^(-iy UoAq) W 3 {q) dq. 



3=0 



(34) 



We will now simplify the second line of this expression for K n+ \. The last sum can be seen as 
an application of Laplace's theorem: 



n— 1 



^{-l) j uojiq) Wj{q) = (-l) n+1 uo,n(q)W(uo,o,-,u ,n-i)(q), (35) 

3=0 



10 



that is, the determinant W(uo,o> ■■■■>UQ,n— l) is expanded with respect to its first row. Next, we 
combine (|35p with its two factors in front, expressing u n ^ n by means of (|25p and ([5]): 



U n ,n(q) 



which leads to the following result: 



h(q) 



2 W(u 0i0 ,...,uo t n)(q) 



(36) 



Kv) 

m 



Un,n(q) 
W(q) 



n-1 



\n+l 



j=0 



After substitution of this expression into (|34p . we get the following form of our propagator 
relation: 



K n+1 (x,y,t) 



Ky) 

m 



n + 1 



n-1 



(-1) £n+i,o2J-i) J -^y 



K (x,q,t) u 0) j(q) dq + 



(v,b) 



'%) 1 



/(y) «n,n(2/) 



L n+1,0 y K o{x,q,t) u , n (q) 



dq. 



(y,b) 



Using once more our representation ([36]) for u n)n , we obtain 

I n-1 



K n+1 (x,y,t) 



Kv) 



(_1V» T x W-lV Wj ^ 



3=0 



K (x,q,t) UQj(q) dq + 



(v,b) 



+ 



L n+i,o / K (x,q,t) u ,n{q) dq. 



f{y)\ W{u , ,--,u ,n)(y) 



The second line of the right-hand side turns out to be the n-th term of the sum in the first line, 
such that we get the following final form of our propagator relation: 



K n+1 (x,y,t) 



n+1 



(-1) L n + i,o2j-iy -^y 



K (x,q,t) u j (q) dq. 



(y,b) 



This expression coincides with (I29p if n is set to n + 1, and the induction is complete. 



3.3 Isospectrality 

We will now consider the remaining case of a SUSY transformation that renders the discrete 
spectrum of the initial boundary- value problem (Q]), © and its final counterpart Q, Q the 
same. In particular, we assume that in each transformation step of the chain the discrete 
spectrum is preserved. Consequently, the auxiliary functions Uq, u n -i used in the SUSY 
transformation ([5]) fulfill only one of the boundary conditions ([2]). Adopting the settings from 
[9], from now on we will require our boundary- value problem to be defined on the whole real 
line, that is, the quantities a, b in ([2]) stand for negative infinity and infinity, respectively. For 
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the sake of simplicity let us first assume that our auxiliary solutions Uq, ...,u n -i only fulfill the 
first boundary condition, that is, 



lim Uj(x) 



0. 



(37) 



The general case of some Uj fulfilling the first boundary condition, and some fulfilling the second 
boundary condition will arise easily once we have established our findings for the setting (|37p . 
More precisely, we will now prove that the following relation between the propagators Kq and 
K n of our boundary- value problems ([!]), (J2j) and ([3]), (UJ) holds: 



K n (x,y,t) 



(-iy 



f(y)\ 



n— 1 y 
j = ® — oo 



(z) cfa. (38) 



The proof of this propagator relation will follow the same steps that were taken in [9] . We must 
show that K n satisfies the generalized time-dependent Schrodinger equation associated to the 
transformed boundary- value problem ([3]), dH), and that for t = the propagator (f38|) becomes 
a delta. In order to prove this latter statement, let us evaluate (f38l) for t = 0. Introducing the 
Heaviside distribution 9, we have 



K n (x,y,0) 



(-iy 
(-iy 
(-iy 
(-iy 



Hy) 

m 

Hy) 
f(y) 

Hy) 
f(y) 

%) 

m 



n-l 

3=0 

n-l 

n-l 



W n (7/) 

W(y) 

Wn(y) 



iCo(^, z, 0) nj(z) 
6{y — x) Uj(x) 



i=o 

n-l 



W(y) 

n 

Wn(y) { f(x)\ 2 W{uQ,...,u j - 1 ,e{y-x) Uj )(x) 



W{y) \h(x 



W(u , ...,uj-i)(x) 



(39) 



Let us now rewrite the Wronskian that involves the Heaviside distribution. Note that the fol- 
lowing argument is the same that was used in [9], we include it here for the sake of completeness. 
In order to rewrite the last Wronskian in (j39|) . we first need to consider the derivatives: 



d" 



dx r - 



6(y — x) Uj(x) 



m—l 

^ C-km 
k=0 



d 



m~k 



dx 



m—k 



: 0(y-x) 



uf\x) + e(y-x)uf l) (x), (40) 



(m) , 



where the Ckm denote constants. If we use this expression to replace the derivatives in the last 
Wronskian of f|39|) and convert the sum in (|40p into a sum of two Wronskians, then the last term 
vanishes, as W(uq, ...,u n -i,Uj) = for j = 0, — 1. We then arrive at 



W(u ,...,Uj-i,6(y - x) Uj)(x) 
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/ u (x) 
u' (x) 



det 



n n _i(x) 

u' n _x(x) -5(x-y)uj{x) 



\ 



( n ) < \ 
u '{x) 



U n-l( X ) £ c kn 
k=0 



n-1 



dx r ' 



9{y — x) ) u*p {x) 



• (41) 



Let us point out that each element of the last column is the sum given by the entry in the lower 
right corner, evaluated at n = 0, 1,2, .... Evaluation of the determinant in (|4"Tj) with respect to 
the last column gives 



W(uo,...,Uj-i,9(y -x)uj)(x) = 

n 

= (-1)" ^T(-l) m W nm (x)Y,C k r 



m—1 



m=0 



k=0 



d 



m—k 



dx m k 



: 0(y - x) 



uf\x), (42) 



where as usual the W nm denote determinants of the minor matrices obtained by deleting the 
n-th row and the m-th column. We will now substitute this result into our propagator (|39p and 
show that it gives a delta, which is equivalent to 



K n (x,y,0) (p(x) dx 



(43) 



for all admissible test functions <j), recall that these functions are smooth and vanish at the 
infinities. After substitution of (j4~2|) into (j39|h we obtain 



K n (x,y,0) <j)(x) dx 



n— 1 n m—1 



W 3 (y) 



j=0 m=l k=0 



W(y) 



d 



dx 



m—k 



: 0(y-x) 



t W nm (x) f(x) uf\x) 
W{x) 



dx. (44) 



The derivatives applied to the Heaviside function can be removed by using 0' = 8 in the distri- 
butional sense. Taking into account 5(y — x) = 5(x — y) and applying the definition of the delta 
function's derivative, we obtain from 



K n (x,y,0) <j)(x) dx 



Hy) 

m 



-, 21 n—1 n m—1 



WAy) 



j=0 m=l k=0 



Qm—k—1 

d^j-i^y-^) 



h{x) 

W) 



W(y) 

§ W nm (x) f{x) uf\x) 



W(x) 



dx 



Hy) 
f(y) 

_&_ 

dy 



n—l n m—1 



WAy) 



m—k—l 



%) 

f(y) 



W(y) 

t W nm {y) f(y) uf\y) 



j=0 m=l k=0 
k \ ( rt/.ATS w t(„.\ ..(*) 



W(y) 



dy. 



(45) 
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We now apply the Leinniz rule to the last term on the right-hand side and make use of the rule 



n-1 
3=0 

which turns (|45p into 

K n (x,y,0) <j>(x) dx 



W(y) U J {V) 



(-1)' 

f(y), 



W nn (y) f(y) 
W(y) 



n-1 



k=0 



(46) 



because W nn = W and the sum in (|4"B]) equals (— l) n , see e.g. [9]. Hence, we have shown that 
K n (x,y,0) = S(x — y). It remains to prove that the propagator K n solves the time-dependent 
Schrodinger equation associated with our generalized boundary- value problem ([3]), (HD, with 
respect to both spatial variables x and y. This is obvious in the first K n is obtained 

from Kq by application of L^ , which maps solutions of our initial boundary- value problem ([JJ, 
([2]) onto its transformed counterpart (|3|), @. Since Kq solves the initial problem, it follows that 
K n must solve the transformed problem. Regarding the second variable y we must substitute 
the explicit form of K n , as given in (j38|) into the time-dependent Schrodinger equation and show 
that it it fulfilled. This short calculation follows exactly the same steps as in [9j, such that we 
omit to show it here. Let us now consider the case where the auxiliary solutions in our SUSY 
transformation fulfill the second boundary condition in ([2]), that is, 



lim Ujix) 

x— >oo 



0. 



In this case, one uses the same argumentation as given above and arrives as the propagator 
relation 



K n (x,y,t) 



-1) 



n-1 



f(y) 



L*n,0 U(-l) 



Kq(x, z, t) Uj(z) dz. 



Finally, if the first M auxiliary solutions fulfill the first boundary condition, and the remaining 
N — M auxiliary solutions satisfy the second boundary condition, then our propagator relation 
will read 



K n (x,y,t) = (-1)" 

+ (-1)"" 1 



M 
M 



TX 



M-l 

n,0 E(" X ) 
3=0 

n-1 

L n,0 E( 
j=M 



Wj(y) 
W(y) 



-IV 



Wj(y) 
W(y) 



Ko(x,z,t) Uj(z) dz + 



Kq(x, z, t) Uj(z) dz. 



As before, it is straightforward to see that our propagator relations reduce correctly to the 
conventional ones constructed in [9], if we set / = h = 1. 
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4 Concluding remarks 



We have shown that the relations between propagators of SUSY-linked Schrodinger equations 
extend to the linearly generalized case. In particular, all possible SUSY scenarios (creation, 
annihilation of spectral values and isospectrality) have been verified and found to match corre- 
sponding findings in [9]. 
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